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Abstract 



It is found here that the l/2 + first excited state of 9 Be is a virtual state 
with the energy of -23.5 KeV. The line shape for the excitation of the state 
is approximated with a simple analytic form based on the effective range 
expansion. The partner in 9 B of this state is found to be a resonance with a 
maximum in the peak at about 1.1 MeV, FWHM of 1.5 MeV, and complex 
energy of 0.6 — i0.75 MeV. The line shape for its excitation is calculated in 
terms of the p— 8 Be phase shift. The phase shifts are obtained from N— 8 Be 
effective potentials deduced from the data on the photodisintegration of 9 Be. 
A possibility for direct extraction of the energy of the resonant state from 
experimental data is also discussed, and an expression for a residue at a virtual 
state pole in terms of a quadrature taken over the virtual state eigenfunction 
is given. 
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I. INTRODUCTION 



In the present work, we obtain the properties of the first excited states of A = 9 nu- 
clei. At low energy these nuclei provide a clean example of three-cluster systems with 
not easily distortable constituents, thus serving as a test ground for theoretical multiclus- 
ter approaches, see e.g. . Reliable information on the properties of these nuclei would 
therefore be very valuable. However, extraction of this information directly from the avail- 
able experimental data is hampered by their ambiguities. The corresponding data on the 
low-energy photodisintegration of 9 Be are not in mutual agreement, see [|J. All states in 
9 B are particle unbound which hinders the search for excited states. There exists a long- 
standing controversy concerning the properties of the 9 B first excited state, as obtained both 
experimentally and theoretically, see @||. Theoretical models can help to analyze the data 
and also guide future experiments. Along these lines, in Ref. || a semi-microscopic model 
to describe the low-energy photodisintegration of 9 Be has been developed. An estimation 
of the reliability of various data sets has been obtained with the help of this model, and a 
theoretical photodisintegration cross section has been derived for astrophysical applications. 
The model provides N— 8 Be effective potentials that reproduce the energy dependence of 
the cross section. Using them the properties of the first excited state of 9 Be are extracted 
below and estimates for its 9 B partner are obtained. 

In the next section we elucidate the nature of the first excited state of 9 Be and obtain 
its position. We obtain an analytic form for the line shape of the state and investigate to 
what degree the line shape is independent of the specific excitation process. In Sec. 3 the 
first excited state of the 9 B nucleus is studied. The expression for the line shape is given, 
and the position, width and the complex energy of the state are calculated. A possibility to 
obtain the latter quantity from a direct fit to the line shape is also discussed. In Sec. 4 our 
results are discussed along with those in the literature. 

Our considerations on the shape of the line in Sec. 2 and 3 are of a rather general 
applicability. In the Appendix a formula expressing a residue at a virtual state pole in 
terms of a quadrature taken over the virtual state eigenfunction is given. 

II. THE FIRST EXCITED LEVEL OF 9 BE 

We proceed from the dynamic input for the description of the system obtained in Ref. [[3] . 
In that work, the 9 Be photodisintegration cross section has been calculated in the framework 
of the following model. The three-body a + a + n representation of the system has been 
adopted. The 9 Be ground state wave function has been calculated from the three-body 
dynamic equation with aa and an potentials. The final continuum state has been chosen as 
a product of the intrinsic wave function of the 8 Be resonance and the n— 8 Be relative motion 
function. The latter wave function has been calculated from the Woods-Saxon potential 
whose parameters were determined by fitting to several radioactive isotope data. 

The model thus gives n— 8 Be effective potentials that allow extracting the properties of 
the first excited state of 9 Be. First, let us comment on the status of the results obtained 
in this way. The model accounts for only the two-body n— 8 Be photodisintegration and 
disregards the direct three-body a + a+n disintegration. The cross section for the direct a + 
a + n disintegration is presumably very small due to the threshold regime which is confirmed 
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in experiment, see 0. Further, the properties of the state considered are determined by the 
phase shift 5 of n— 8 Be scattering. So, the question arises whether this phase shift can be 
correctly obtained with the help of the model. In this connection one can admit that in the 
vicinity of the first excited state of 9 Be the energy dependence of the photodisintegration 
cross section is also determined by this phase shift. To a degree this holds true, use of 
the two-body 8 Be+n dynamics to extract the properties of the state does not contain any 
restrictions. In Ref. j3| the cross section has been fitted for excitation energies up to 0.5 
MeV, and in this region the corresponding energy dependence sm 2 S/k dominates the cross 
section, see Fig. 1 below. Moreover, the initial 9 Be ground state is described realistically in 
this model. Therefore, use of the two-body dynamics seems to be sufficient for extracting 
the properties of the state considered. In general, these properties can also be extracted via 
the fit of an assumed line shape directly to the data. However, use of the two-body dynamics 
in conjunction with shell-model considerations allowed selection between alternative data 
sets in Ref. ||. Below it will also help us to obtain the properties of the analog 9 B state 
relying on the 9 Be photodisintegration data. 

The fit to the photodisintegration cross section || provides the potentials 

V(r) = Vb/[lW r -* )/6 ] (1) 

describing n— 8 Be scattering. The best version 

V = 35.99 MeV, R = 3.126 fm, b = 0.8108 fm (2) 

will be used below. Other potential versions lead to similar results as commented below. 

Let us denote E = (hk) 2 /(2fi) the energy of the relative n— 8 Be motion. The potential 
(H), (||) leads to a n— 8 Be scattering amplitude f(k) having a pole at k — —in, k > 0: 

'<*> - irk (3) 

as k — > —in. This means that the first excited state of 9 Be is a virtual state strongly coupled 
to the n— 8 Be channel. The energy of the state is 

- (hK) 2 /(2fi) = —E = -23.53 KeV. (4) 

This value was obtained via solving the eigenvalue problem with the virtual state boundary 
condition, cf. the Appendix. 

Let us consider the shape of the line for a transition proceeding via the virtual state. 
The shape depends partly on the specific way of excitation of the state. The main energy 
dependence is universal, however, and it is determined by the intrinsic properties of the 
state. We shall compare the universal energy dependence to the energy dependence for the 
process of 9 Be photodisintegration and thus estimate the dependence of the line shape on 
the specific process. 

We assume that the cross section for formation of 9 Be in the continuum with the energy 
E in the vicinity of the excited level can approximately be presented in the form 

a ~ J k 2 dk\(^- k (J = l/2 + )>| 2 <f(fc 2 /(2/x) - E). (5) 
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Here $ is some localized state in the subspace of the 9 Be degrees of freedom. It includes the 
g.s. of 9 Be and the (effective) transition operator accounting for the excitation mechanism 
(photodisintegration, e.g.). The transition operator may have an energy dependence that is 
smooth at E — > 0. The function vpjjT is the continuum spectrum function that corresponds to 
outgoing n and 8 Be fragments in the relative s-state thus describing the two-body 8 Be+n 
photodisintegration. The representation (§) implies that the formation process can be de- 
scribed in the subspace of the 9 Be degrees of freedom. Under this condition properties of a 
state reveal themselves independently of the interactions with other particles participating 
in its formation.^ In Eq. (|5|) the normalization condition \^k') = — k')/k 2 should be 
fulfilled, so that the n — 8 Be relative motion function in vp^T at large distances is normalized 
to e~ lS sin(fcr + 5)/ (kr). 

We may write the fc-integral of (5) as a contour integral where the pole contribution is 
given by (3), neglecting influence of other degrees of freedom in <3>, ^ . 

The energy dependence of the cross section (|5]) in the vicinity of the virtual level is 
sin 2 6/ k. This energy dependence arises if one replaces the outer part sin(A;r + 5)/(kr) of 
the n— 8 Be relative motion function in the matrix element from Eq. (|^) by sin 5/ (kr). This 
can be done under the conditions (kR) 2 <C 1, and R <C \a\. Here a is the scattering length, 
and R is chosen such that the relative contribution of r > R values to Eq. (|5|) is small. 
R exceeds the range of the n— 8 Be interaction, and the inner part of the final state wave 
function acquires the same energy dependence as the outer one. The sin 2 S/k energy 
dependence is nothing else as the so called Migdal-Watson factor ||. 

We can rewrite the corresponding contribution to the cross section as const-/c|/(A;)| 2 
where f(k) is the s-wave n— 8 Be scattering amplitude. Almost the same accuracy is kept if 
one takes f(k) within the effective range approximation 

f( k ) = [-l/a + (l/2)r k 2 -ik}- 1 . (6) 

Then one obtains \f(k)\ 2 ~ [(E + E)(E + Ex)]' 1 . In the model ©, (§) E x = 1.569 MeV. 
Using this expression, it is convenient to present the cross section in the following form 



fiA 2VEEE + E, 

Since E E\ the energy dependence of the cross section from the threshold up to its 
maximum is entirely determined by the second factor in Eq. (|?|). It is then seen that the 
maximum of the cross section occurs practically at the E value. With a high precision a m 
is the value of the cross section at the maximum. As one can see from Eq. (||) the lowest 
order correction to Eq. (^]) is of the form 1 — (E/E 2 ) where E 2 is not small and depends on a 
specific excitation process. To a certain degree it can be accounted for via a renormalization 



1 In the framework of the model Q only the two-body 8 Be+n channel is retained in the function 
while the incoming three-body a + a + n channels are disregarded. This model assumption 
is also inherent to all the previous work. One can see, however, that to a first approximation this 
does not change the energy dependence of the cross section and thus does not influence the results. 
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of E±. Eq. (HD is the Breit-Wigner formula with T = T(E), (or, the one-level R matrix 
expression) rewritten in a different form. 

In Fig. 1 various universal expressions for the shape of the line are compared with the 
exact line shape for the photodisintegration process calculated in the model of Ref. [[|. 
(E = E 1 — E t h-) The full curve represents the latter line shape, the long-dashed curve is 
the energy dependence ~ k\f(k)\ 2 , the dash-dotted curve represents the expression fl7|), and 
the dotted curve is the energy dependence given by the second factor from Eq. (H). The 
FWHM value for the photodisintegration process provided by the first mentioned curve is 
196 KeV. The FWHM values provided by the next two curves are 230 and 240 KeV. The 
latter two values are process-independent. Thus FWHM may depend rather sizably on the 
excitation process. On the contrary, the position of the maximum in the peak is practically 
process-independent coinciding with the absolute value of the energy of the virtual state. 

Experimentally, the energy dependence E X I 2 {E + E)" 1 , specific to a virtual state could 
be confirmed by measuring the shape of the line from the threshold up to the region of the 
maximum with the tagged photon techniques. A simpler, while indirect, way is to study the 
whole peak in the (e, e') reaction and extract the properties of the state in the way similar 
to that used above for photodisintegration. In case of accurate and detailed (e, e') or (p,p') 
data the energy E of the state could also be extracted from the fit of the Eq. (pi) form or 
its above-mentioned extension. 

Besides the energy E, or the pole position of the n— 8 Be scattering amplitude, another 
quantity to be reproduced in a microscopic calculation is the residue Co in the pole, Eq. @. 
The exact Co value was calculated with the formula derived in the Appendix that represents 
Co as an integral taken over the virtual state eigenfunction. The result is cq = 0.7837. We 
note that with high accuracy both E (or k, see (||),) and c can be expressed in terms of the 
n— 8 Be scattering length a and the effective range r . From Eq. ([]) one obtains 

K = r - 1 [(l-2r /a) 1 / 2 -l], 
c = (l-2r /a)- 1/2 = (1 + ^ro)- 1 . (8) 

The effective range parameters for the potential ([]]), @ are 

a = -27.65 fm, r = 8.788 fm, (9) 

and the expected accuracy of the above expressions is about |ro/a| 3 ~ 1%. Their numerical 
values proved to be even more accurate: E = 23.51 KeV, Co = 0.7819. As it is explained 
in [U all other acceptable potentials found there lead to the a and r values very close to 
those in Eq. (|). Hence the properties of the virtual state given by these potentials are 
quite similar. 

III. THE FIRST EXCITED LEVEL OF 9 B 

To calculate the resonant peak for an excitation of the 9 B l/2 + level, we proceed again 
from Eq. (|5]), only with a different expression for the pole term but still assuming that there 
is only one decay channel, p— 8 Be, for the state considered. In the 9 B case, disregarding 
the possibility for the decay into the N + a + a channel, being less substantiated than in 
the 9 Be case because of a higher energy with respect to the three-body threshold, seems 
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still to be reasonable. The n— 8 Be potentials obtained in Ref. || lead presumably to a good 
description of the n— 8 Be l/2 + phase shifts. Then one may suggest that the p— 8 Be l/2 + 
phase shifts will also be properly reproduced with these potentials with an addition of the 
Coulomb interaction. This will suffice to obtain approximately the main properties of the 
9 B 1/2+ level. 

Beyond the range of the nuclear potential the p— 8 Be relative motion function x(r)//cr 
entering ^ in Eq. (|5|) may be represented in the form 

X = Gsind + Fcosd = (sin 5 / C)[C (G + F cot 5)] = (sin 8/C)<j)(k, r). (10) 

Here 5 is the nuclear phase shift, C 2 = 27n][exp(27n]) — 1] _1 is the Coulomb penetrability 
factor, and F and G are the Coulomb functions. The possible resonant pole is contained in 
the factor sin 5 while the function defined in Eq. QIUD is smooth at low energies. (At k — > 
4>(k, r) — > zKx(z) — (a c /Aa)zIx(z), z = 2(2r/a c ) 1//2 , where a is the scattering length, and a c is 
the Bohr radius.) To a first approximation one may use 4>(ko, r) in the calculation of the cross 
section (|5p, fco being chosen in the vicinity of the resonance. Hence the energy dependence 
of the excitation cross section in the resonant peak region is given by the expression 

sin 2 6 _ \f(k)\ 2 _ kC 2 
kC 2 C 2 ~ [kC 2 cot 5) 2 + {kC 2 ) 2 

Here / is the nuclear amplitude of p— 8 Be scattering. Thus, quite naturally, up to process- 
dependent corrections the resonant cross section for the excitation of the state is propor- 
tional to the scattering cross section times a universal factor proportional to a width. The 
quantity kC 2 cotS in (|TTD allows the well-known representation [— 1/a + (l/2)ro/c 2 + ...] — 
(2/a c )h(l/a c k), showing that it is smooth and finite at E — 0. The kC 2 behavior of the 
cross section at E — > is seen directly from Eq. ([5]) and the properties of the Coulomb 
functions. 

The p— 8 Be phase shifts entering Eq. ( [XT]) were obtained from the Schrodinger equation 
with the nuclear potential ([!]), @ plus the Coulomb interaction. The latter took into account 
the density distribution in 8 Be and the charge distribution in the a particles: 

V C oui(r) = (Ae 2 /r)(2/7T) (sin qrj q)F a (q)I(q)dq, (12) 

Jo 

where F a (q) is the charge form factor of the a particle 0, and 

I(q) = / (2/qp)sm(qp/2)^ 2 (p)dp, 
Jo 

ip being the 8 Be wave function. 

The energy dependence flTT|) obtained indeed proved to exhibit a pronounced peak, and 
it is shown in Fig. 2 with a solid curve. The position E max of the maximum in the peak 
with respect to the p— 8 Be threshold and the FWHM value are the following: E max =1.13 
MeV, FWHM=1.64 MeV. The long-dashed curve in Fig. 2 represents the peak for the other 
acceptable nuclear potential of the Woods-Saxon form found in Ref. whose parameters 
are 

V Q = 52.86 MeV, R = 2.006 fm, b = 1.051 fm. (13) 
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(11) 



For this case E max = 1.02 MeV, and FWHM=1.43 MeV. 

If the peak obtained corresponds to a state then the scattering amplitude should have a 
pole in the vicinity of E max . Thus we shall search for such a pole. One can show that when 
the long-range Coulomb interaction is switched on, virtual states cease to exist, turning 
normally to complex-energy resonances, and this applies to our case. (It applies also to, 

e.g., p p scattering where, in contrast to what is often said, there are no virtual states.) 

The energy E = Eq — iT /2 of a resonance being connected to the pole of a scattering 
amplitude, is process-independent and thus characterizes a resonance quite precisely even 
in case of a broad width. 

In our case this quantity is computed as a complex eigenvalue in the p— 8 Be Schrddinger 
equation. Here we used the codes of Ref. || or, when this is inoperative, the p-series 
(14.1.4) for F and the corresponding expansion (14.1.14) - (14.1.19) for G from Ref. ||.[] 
These expansions are fast convergent for not extremely high values of \p\, \rj\. 

The resonant pole was found with the following parameters: Eq = 0.60 MeV, T/2 = 0.77 
MeV, and E = 0.56 MeV, T/2 = 0.70 MeV for the potentials (§) and flT3D, respectively. 
The positions Eq of the resonance are shifted downwards with respect to the maxima in 
the peak found above, while the widths T are close to the FWHM values. Similar trends 
were observed for 5 He, 5 Li resonances pD[. It is interesting to note that the height of the 



Coulomb barrier in our case proved to be 0.63 MeV only, i.e. the resonance with a finite 
width sits at the top of the barrier. This is possible only when the width of the resonance 
becomes so broad that it is comparable to Eq. 

In conclusion, let us add the following two comments. First, since the energy E = 
Eq — iT/2 of the resonance is a convenient process-independent quantity to compare with 
theory, it would be very useful to extract it directly from future experimental data without 
constructing a model for the excitation process. Such a task has been accomplished suc- 



cessfully in some cases, see e.g. pO| , pT| . However it is not clear whether this is possible in 
practice for the broad resonance we consider. We perform a numerical experiment to clarify 
this issue. We explore the possibility to reconstruct the Eq and T values obtained proceeding 
from a reasonable fit to the resonant peak of the form (|IID , see Fig. 2, calculated in the 
same model, Eq. (0). We fit the peak with the expression 

Gl [a 2 + a 3 (E - E>) + a 4 (E - E') 2 ] 2 + (kC 2 ) 2 (M) 

related to Eq. Oil]). Here E' = 1.5 MeV, the fit is extended over the values (E' — T/2) < 
E < (E' + T/2) at which the cross section exceeds a half of its value at the maximum, and eij 
are fitting parameters. At choosing the form of the expression flli]) it was taken into account 



2 We rewrite the latter expansion for / = in the form valid for complex rj values: 

oo 

G( V ,p) = C^{2rip^(ri,p)[ln(2p) + (l/2)^(i V ) +^(-ir,)} + 2 7 - 1] + Y l a k(v)p k }- 

k=0 

Here ip(z) = T'(z)/T(z), 7 is the Euler constant, and the quantities $(77, p) and 0^(77) are defined 
in [@. 
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that the effective range approximation is not accurate enough in the peak region and that 
kC 2 cot 5 has no minimum in that region in our case. The resonant E value, E = E — iT/2, 
is obtained as a zero of the expression a 2 + a 3 (E — E') + a 4 (E — E') 2 — ikC 2 that corresponds 
to the pole of the scattering amplitude. 

The search of the least square minimum was performed with respect to Eq, T, and as 
with values taken on some grids. At given values of these parameters a 2 and were fixed 
via equating to zero real and imaginary parts of the above expression. The a\ value was 
found analytically from the least square minimum requirement. (If kC 2 cot 5 is completely 
reproduced by the fit, we should have a\ = 1.) 

It occurred that the fit is unstable in this form, and a\ takes unrealistic values. Let us 
then remove the search of the T value admitting instead the hypothesis T =FWHM. This 
condition is fulfilled with sufficient accuracy in our case, and also in the 5 He, 5 Li cases JTU . 
This modified procedure leads to E values that are quite stable and close to the true one. 
The stability was checked by changing the stepsize in the a 3 parameter. 

The second comment concerns the value of the effective range r for the p— 8 Be scattering. 
For the potential (0), for example, it proved to be 3.01 fm i.e. much smaller then that from 
Eq. (H) for the n— 8 Be case. This, however, does not mean that the range of the p— 8 Be 
nuclear interaction is sizably different from that for the n— 8 Be interaction. The effective 
range can serve as a measure of the interaction range when a ^> r and, in addition, the 
zero energy scattering wave function u(r) entering the effective range definition is close to 
unity at the edge of the well R. This holds true for the neutral case but not when the rather 
strong Coulomb interaction is present. The latter interaction suppresses the wave function 
at the R value, so that the range of nuclear interaction may be estimated as ro/u 2 (R). 



IV. DISCUSSION 

Basing on the model of Ref. || we obtained that the first excited level of 9 Be is a virtual 
state with the energy — E =-23.5 KeV with respect to the n— 9 Be threshold. It is shown 
that the peak position for an excitation of this state practically coincides with the E value 
while the FWHM value of the peak may sizably depend on a specific excitation process. 
Within the effective range approximation the properties of the state are reproduced with 
a high accuracy. The line shape is aproximated with a simple analytic form based on the 
effective range expansion. 

In Ref. a single-level .R-matrix fit to the photodisintegration data of Ref. [13] led 



to the description of the state considered as a complex-energy resonant state. However the 



data [13J are at variance with the earlier data, and it was concluded in Ref. [y] that the 
latter ones are preferable. The positions of maxima in the available experimental 9 Be(p,p') 
spectra ]|14| , |T5|] do not contradict the above listed value of 23.5 KeV. 



In Ref. JT5[ the quantities pertaining to the l/2 + state of 9 Be were presented in the form 
of parameters entering the Breit-Wigner formula 

b(ei e ) T JMilL 

9 he 7 K ' i } {E-E R Y + \r(E)/2Y< 



<E) = — -E,B(E1, E,) /rt F ; 2 " wo 12 » ( 15 ) 



that was fitted to the (e, e') data of that paper. Here T(E) = GyE, G being the reduced 
width of the level. The values Er — 19±7 KeV, V = 217±10 KeV were reported, the latter 
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value refers presumably to T(Er). These values are quoted in the review article |16[| . They 



cannot be correct since they lead to a quite unrealistic E value, i.e. that of the maximum 
of an excitation cross section, of 0.6 KeV. Based on Fig. 6 from Ref. IT3] one may suggest, 
however, that these values do not refer to E R and Y from Eq. (|i~5"D but to the position of the 
maximum of the cross section and the FWHM, respectively. If it is the case, these values are 
in agreement with ours. (Absence of the information on energy dependence of B(E1,E 7 ) 
inhibits a precise conclusion concerning the FWHM.) 

In Ref. a microscopic study of the spectrum of 9 Be at the 9-nucleon level was un- 
dertaken. The l/2 + level was not detected at all. The reason may lie in that the method 
of complex scaling used in that work is suited for search of complex-energy resonant states 
but not virtual states. 

Further, we studied the first excited state in 9 B in the present work using the nuclear 
potentials derived in || from the analysis of the 9 Be photodisintegration data and adding the 
Coulomb interaction between p and 8 Be. An approximate process-independent expression 
for the line shape of the state in terms of the p— 8 Be phase shift is obtained. The position of 
the peak with respect to the p— 8 Be threshold and its FWHM given by this expression are 
about 1.1 MeV and 1.5 MeV, respectively. The peak considered is caused by a complex- 
energy pole in the p— 8 Be scattering amplitude. The pole proved to be located at energy 
E = E — iT/2 with E ~ 0.6 MeV and T ~ 1.5 MeV. 

In Ref. ||] a prediction for the peak of the 9 B resonance was obtained: E max = 1.13 MeV, 
FWHM=1.40 MeV. It is close to ours while the underlying assumptions of that work were 
rather different. Our common features with that work consist in use of p— 8 Be dynamics to 
obtain the resonance and in obtaining the p— 8 Be Woods-Saxon potential from 9 Be photo- 
disintegration data. The big differences consist in the conditions from which parameters of 
the potential are deduced, in the parameters themselves, and in a prescription to calculate 
the resonance. In Ref. || the two-body dynamics were used both for the inital ground state 
and the final continuum state of 9 Be. Depths of the potentials were varied whereas their 
range and diffuseness were kept at their "classic" vaues. The data of Ref. |L3| were fitted 



(probably up to a normalization). The fit was only moderately good. The line shape of the 
9 B resonance was computed as if it were excited due to a fictitious dipole transition from 
the ground state of 9 Be. Our potentials were deduced || at the assumption of three-body 
a + a + n dynamics for the ground state of 9 Be and two-body dynamics for its continuum. 
Range and diffuseness of the potential were varied in addition to the depth. The data of 
Ref. were concluded to be less preferable, and the alternative earlier data were fitted 
without freedom in an absolute normalization. The fit is statistically quite good. The line 
shape of the 9 B resonance was computed from Eq. (|TT|). The Coulomb interaction was 
treated more accurately, and also the pole position E — iT/2 was calculated. Keeping range 
and diffuseness of the potential at their " classic" values was perhaps an important point in 
the analysis of Ref. |4j] which allowed obtaining correct results even at assumptions that are 
not completely valid. An interesting point is that the different prescriptions for calculating 
a resonance in our work and in Ref. lead to similar results. To check this we calculated 
the resonance for the potential of Ref. [[| with the prescription of Eq. (0). The results 
are E max = 1.06 MeV, FWHM=1.47 MeV, and they are close to those reported in Ref. J§. 
(The difference in treatment of the Coulomb interaction should also be taken into account 
here.) 
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In Ref. []17] the energy of the l/2 + first excited state of 9 B was derived from i?-matrix 
parameters fitted [O] to the 9 Be data of Ref. [13] and from values of the Coulomb displace- 



ment energy calculated with the help of the shell model. The result E ~ 2 MeV differs 
considerably from that of our work and Ref. 0] leading to the inverted value of the so called 
Thomas-Ehrmann shift. As it is mentioned above in connection with our comparisons to 
Refs. fl"2]]4jj, use of this particular set of data to derive i?-matrix parameters might be a 



disadvantage. 

We are indebted to M.V. Zhukov for useful discussion in the course of this study. A 
part of this work was done during the stay of V.D.E. at the Niels Bohr Institute, and he 
expresses his gratitude for the kind hospitality. The work was partially supported by Russian 
Foundation for Basic Research (grants 96-15-96548 and 97-02-17003). 

APPENDIX: RESIDUE PROPERTIES OF VIRTUAL STATE POLES 
OF A SCATTERING AMPLITUDE 

In the vicinity of the virtual state pole the two-body scattering amplitude behaves ac- 
cording to Eq. (|). We obtain a formula that expresses a residue Co in terms of a quadrature 
taken over the virtual state eigenfunction. The result is c = (—l) l /I v , where I v is given 
by Eq. ( |A4| ) below with A(r) entering the virtual state eigenfunction ( |A3| ) . The latter 
function is readily calculated from the Schrodinger equation. (The consideration is formally 
applicable to a motion in a central potential with any orbital momentum I although the 
/ > virtual state case is of little interest.) 

Let ip{r) be the solution to the Schrodinger equation regular at the origin. Be- 
yond the range of a potential the function x(r) = rtf)(r) behaves as exp(ikr) + 
(—l) l+1 S~ 1 (k)exp(—ikr). Complex k values are allowed. A virtual state corresponds to 
a pole of the S'-matrix S(k) at k = —in, k being real and positive.^ In terms of n, 
x{r) = exp(«;r) + (— 1)' +1 S _1 exp(— nr) at large r, so the boundary condition for virtual 
states consists in absence of the decreasing exponential. If the eigenvalue is sufficiently 
small then the eigensolution can be found directly using a logarithmic derivative at suffi- 
ciently large r. This holds true in our case. 

For bound states, the well-known formula corresponding to that we want to derive is 



the following |TR[T9]. Let the bound state wave function be normalized to unity, and A the 



coefficient in its asymptotics Aexp(-nr). Then 



as k tends to in. 

In the virtual state case we shall proceed from the relation 



dx dx d 2 x 2m 
dr dE ^ drdE 



5/ X\r')dr>. (A2) 
n Jo 



Tt is implied that at large r the potential decreases more rapidly than exp(— nr), cf 
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It is similar to that used in Ref. || for the derivation of Eq. ( |A1|) . We take r to be 
large enough so that in the vicinity of the virtual state pole x{ r i E) — exp(Kr) + a(E + 
E) exp(-Kr). Here E = (hn) 2 /2m. The pole term in the scattering amplitude is then 
(_ iy+ 1 (2,Ka)~ 1 (E + E)" 1 , and the constant a is to be found. We denote Xvij) the virtual 
state eigenfunction, so that Xv = x(E = —E). We set 

Xv(r') = exp(nr') + A(r'). (A3) 



Here A(r) is a rapidly decreasing function. We substitute Eq. ( |A3| ) into the right-hand side 
of Eq. (JA^) and equate the terms of the orders of exp(nr) and unity in both sides of Eq. 
(|A2|) at E tending to — E. Then we obtain 2na = —(2m/h 2 )(I v /2K), with 



roc 

I v = 1 - 2k / [A 2 (r) + 2A(r) exp(/tr)]dr. (A4) 
Jo 

The product A(r) exp(«;r) is a rapidly decreasing function, so the integral converges. Finally, 

/(*) - (AS) 

as k tends to —in, k > 0. 

To display an analogy between the bound and virtual state pole cases, we note that the 
right-hand side of Eq. (|A1|) can be rewritten as 



•IV 



2nNh k — in' 



(A6) 



where Nj, is the norm of the eigenfunction Xb behaving as exp(— nr) at r tending to infinity. 
If one sets Xb( r ) — exp(-nr) + A(r) then 2kN^ takes just the form of Eq. (|A4|) with the 
replacement k —>■ — n. 
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FIGURES 



FIG. 1. Line shapes for transitions to the first excited state of 9 Be. The solid curve represents 
the exact line shape for the photodisintegration process calculated with the model of Ref. [5]. The 
long-dashed curve is the energy dependence k\f(k)\ 2 , where / is the n— 8 Be scattering amplitude. 
The dash-dotted curve stands for Eq. (0). The dotted curve represents the energy dependence 
given by the second factor from Eq. ([?]). 

FIG. 2. Line shapes for transitions to the first excited state of 9 B calculated according to Eq. 
(|Tl|). The solid and long-dashed curves are for the potentials with the parameters @ and (|l3|), 
respectively. 
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